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ABSTRACT 

The LIM homeobox 2 transcription factor Lhx2 is 
known to control crucial aspects of neural develop- 
ment in various species. However, its function in 
human neural development is still elusive. Here, we 
demonstrate that LHX2 plays a critical role in human 
neural differentiation, using human embryonic stem 
cells (hESCs) as a model. In hESC-derived neural 
progenitors (hESC-NPs), LHX2 was found to be ex- 
pressed before PAX6, and co-expressed with early 
neural markers. Conditional ectopic expression of 
LHX2 promoted neural differentiation, whereas dis- 
ruption of LHX2 expression in hESCs significantly 
impaired neural differentiation. Furthermore, we 
have demonstrated that LHX2 regulates neural dif- 
ferentiation at two levels: first, it promotes expres- 
sion of PAX6 by binding to its active enhancers, and 
second, it attenuates BMP and WNT signaling by 
promoting expression of the BMP and WNT antag- 
onist Cerberus 1 gene (CER1), to inhibit non-neural 
differentiation. These findings indicate that LHX2 
regulates the transcription of downstream intrinsic 
and extrinsic molecules that are essential for early 
neural differentiation in human. 



INTRODUCTION 

The neuroepithelium is derived from epiblasts and is the 
origin of the nervous system. Its formation is based on 
the interplay between intrinsic transcription factors and 
extrinsic signaling (1). Experiments in Xenopus and 
rodents suggest that extrinsic signaling molecules, such 



as fibroblast growth factor (FGF), support neural induc- 
tion (2), whereas WNT, SHH, Notch, transforming 
growth factor beta (TGF(3) and BMP hinder neural dif- 
ferentiation (3). On the other hand, intrinsic transcription 
factors, such as Soxl, have been found to control neural 
determination in mice (4). Soxl was further demonstrated 
to influence mouse neural specification by regulating ex- 
pression of neural genes, including Pax6 (5). These 
findings suggest that the differentiation of pluripotent 
stem cells into neuroepithelia is regulated by the interplay 
between signaling factors and neural transcription factors. 

Complicated ethical issues, and lack of a cellular system 
suitable for investigating the molecular mechanisms 
underlying early lineage differentiation, have hampered 
studies into early human development; however, the estab- 
lishment of human embryonic stem cells (hESCs) has 
helped to mitigate these challenges (6). Human ESCs are 
pluripotent, capable of unlimited proliferation in vitro. 
They thus provide an alluring and powerful pipeline 
through which to study the molecular basis of the transi- 
tion from pluripotency to the neural lineage. To date, 
obligated primitive stage neural cells (7,8) and functional 
transplantable neurons have been generated from hESCs 
in vitro under chemically defined neural induction culture 
conditions (9,10). Using this neural induction procedure, 
the genetic factors that contribute to the regulation of 
early human neural differentiation have begun to be 
unearthed (4,5,7,11). 

Transcriptional activator LIM homeobox 2 (Lhx2) has 
roles in neural induction and morphogenesis, as well as 
diverse functions among species (12). Lhx2 homologs are 
found in Drosophila (apterous), zebrafish (bel), Xenopus 
and mouse (12-15), where their sequences are largely 
conserved, but show some divergence in spatial expression 
in the central nervous system (12,16). During mouse 
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embryogenesis, Lhx2 expression initiates before 
neurulation, and it is highly expressed in the neuroepithe- 
lium and ventricular zone (12,17,18). In Lhx2 mutant 
mice, central nervous system defects are detected in the 
dorsal telencephalon, hippocampus, thalamus, optic 
vesicle and olfactory bulb, suggesting that Lhx2 plays a 
critical role in brain morphogenesis (18-26). Although 
recent studies indicate that Lhx2 participates in the deter- 
mination of regional fate decisions specifying cortical 
identity (18,25), the functional role of LHX2 in transcrip- 
tional control of early neural differentiation in human is 
not well understood. 

The expression pattern and role of LHX2 in human 
neural differentiation have not yet been investigated. 
Here, we report that LHX2 is one of the early genes 
switched on in differentiating hESCs, before PAX6 and 
SO XI, and its expression is maintained with that of 
other neural genes in hESC-derived neural progenitors 
(hESC-NPs). Additionally, gain- and loss-of-function ex- 
periments revealed that overexpression of LHX2 
enhances neural differentiation, whereas its knockdown 
impairs neural differentiation in hESCs. We have 
further demonstrated that LHX2 directly regulates ex- 
pression of the neural transcription factor PAX6, and 
the WNT and BMP signaling antagonist Cerberus 1 
(CER1). Taken together, our results demonstrate that 
LHX2 regulates the transcription of downstream intrin- 
sic and extrinsic molecules that control early neural 
differentiation. 



MATERIALS AND METHODS 

Cultivation and neural differentiation of hESCs 

Human ESCs H9 (WiCell, Madison, WI) (6) and NTU1 
(from C.F. Chen, NTUH) (27), and human induced 
pluripotent stem cells (iPSCs) iCFB46 (28), were used at 
a passage number below 50. Undifferentiated hESCs were 
routinely maintained and subcultivated as previously 
described (29). For neural differentiation (Figure 1A), 
hESCs were dissociated with 1 mg/ml dispase (Sigma), 
and then left to form embryoid bodies (EBs) in ultra-low 
attachment dishes (Corning) for 4 days in hESC media 
[Dulbecco's modified Eagle's medium (DMEM)/F12, 
20% knockout serum replacement, 1% non-essential 
amino acids, 2mM L-glutamine, 100 uM P- 
mercaptoethanol and 5ng/ml basic fibroblast growth 
factor (bFGF)]. EBs were transferred for further neural 
differentiation in N2 media [DMEM/F12, N2 supplement, 
1% non-essential amino acids, 2mM L-glutamine and 
20ng/ml bFGF (all purchased from Invitrogen)]. EBs 
were plated on tissue culture dishes (Nunc) in N2 media 
for 6 days of in vitro differentiation (IVD). For gain-of- 
function experiments, cells were separated from feeder 
cells using 1 mg/ml dispase (Sigma) and re-plated on 
Matrigel coated dishes. One day after plating, cells were 
treated with 2 ug/ml doxycyline (Sigma) in DMEM, 3% 
FBS, 1% non-essential amino acid (NEAA) and 2mM l- 
glutamine (Invitrogen) for 3 days. 



Generation of inducible overexpression and constitutive 
knockdown hESCs 

For generation of Tet-On inducible overexpression 
hESCs, EFla and TetR (from pcDNA6-TR, Invitrogen) 
were cloned into pcDNA3.1 as pTetR; and EFla, TetO, 
LHX2 coding sequences cloned from hESC-derived neural 
progenitors, and puroTK (from Y.T. Chen, NTU) were 
cloned into pcDNA3.1 as pTetOLHX2. Electroporation 
of DNA constructs into hESCs and establishment of 
stable clones was performed as described previously (30). 
Briefly, 10 7 undifferentiated hESCs were prepared for each 
electroporation. After pretreatment with Y-27632 
(688000, Merk) for 2h, cells were trypsinized into single 
cells, and used for electroporation or lentiviral infection. 
Electroporation involved 25 ug plasmid and the following 
conditions: 230 V, 500uF. For generation of conditional 
overexpression hESCs, TetR-H9 was established by trans- 
fecting pTetR into H9, followed by drug selection using 
G418 (50 ug/ml; Invitrogen). After establishment of TetR- 
H9, pTetOLHX2 or pTetOGFP was transfected into 
TetR-H9 by electroporation and selected by puromycin 
(400 ug/ml; Sigma) for establishment of iLHX2 or iGFP 
stable clones, respectively. Two iLHX2 clones with higher 
expression levels of LHX2 after 2 ug/ml doxycycline 
(Sigma) induction were selected for subsequent 
overexpression experiments. For generation of constitu- 
tive knockdown hESCs carrying shLHX2, shPAX6, 
shCERl or shLuc, shRNA on pLKO_TRC001, constructs 
were obtained from the RNAi Core, National Science 
Council, Taiwan, and lentiviral particles were produced 
using 293T cells following standard procedures. The 
shRNA constructs were named according to the targeted 
gene and the position of the first shRNA targeting nucleo- 
tide in the gene coding sequence (number). The targeting 
sequences were as follows: shLHX2-831 (GCTTCGGAC 
CATGAAGTCTTA), shLHX2-968 (GCAACCTCTTAC 
GGCAGGAAA), shPAX6-284 (CGTCCATCTTTGCTT 
GGGAAA), shCE7?7-338 (CCAGCCGATAGATGGAA 
TGAA) and shCERl-460 (GCCATGAAGTACATTGG 
GAGA). Lentiviral particles were generated carrying 
shLHX2, shPAX6, shCERl or shLuc. Titers of resulting 
lentiviral particles were examined following standard 
MTT procedures, and 0.01 R.I.U. viral particles were 
introduced into undifferentiated hESCs; cells were then 
selected for with puromycin (400 ug/ml; Sigma) for 
stable clone establishment. After selection, two to three 
single colonies were selected in mixed culture as constitu- 
tive knockdown hESCs. Each stable clone was examined 
for pluripotency gene expression by reverse transcriptase- 
polymerase chain reaction (RT-PCR) and AP staining 
(Supplementary Figure S2). 

Transient overexpression in hESCs 

For transient overexpression of LHX2 and PAX6 in 
hESCs, GFP, LHX2 and PAX6 were cloned into FUW 
to generate FUW-GFP, FUW-LHX2 and FUW-PAX6, 
respectively. Transient transfection of hESCs with 2ug 
FUW-GFP, FUW-LHX2 or FUW-PAX6 was performed 
with Lipofectamine 2000 (Invitrogen), following the 
standard Lipofectamine 2000 transfection protocol. 
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Figure 1. LHX2 is expressed in neural lineage cells derived from hESCs. (A) Procedure of in vitro neural induction from hESCs. (B) RT-PCR and 
(C) qRT-PCR analysis of neural genes using mRNA isolated from differentiating H9 hESCs at the indicated days. LHX2 transcripts could be 
detected earlier (at day 4) than PAX6 and SOX1 transcripts. Error bars represent the mean ± SEM (n = 3). (D) Western blot of proteins isolated 
from differentiating hESCs at the indicated days with specific antibodies against LHX2 and PAX6. (E) Immunocytochemical analysis of 
differentiating hESCs at the indicated days with specific antibodies against PAX6 and LHX2. LHX2 was co-localized with PAX6 in differentiating 
hESCs at various differentiation stages. 
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Reverse transcription-polymerase chain reaction 

RNA was extracted with TRIzol reagent following the 
standard extraction protocol (Molecular Research 
Center). Extracted RNA was reverse-transcribed into 
cDNA with Superscript III Reverse Transcriptase 
(Invitrogen). Each PCR used 25 ng of cDNA. GoTaq 
Green Master Mix (Promega) was used for RT-PCR. 
For quantitative PCR (qPCR), SYBR® FAST 2X qRT- 
PCR Master Mix (KAPA) and a 7900HT Fast Real-Time 
PCR System (Applied Biosystems) were used. Primers are 
listed in Supplementary Table SI. 

Teratoma generation and ethics statement 

For teratoma generation, 5 x 10 6 undifferentiated hESCs 
were subcutaneously injected into NOD-SCID mice. After 
7 weeks, mice were sacrificed and teratomas were isolated. 
All the animal experiments were approved by the Animal 
Care and Use Committee of Academia Sinica, and per- 
formed in accordance with the Institutional Animal Care 
and Use Committee guidelines of Academia Sinica 

Immunocytochemistry, immunohistochemistry and 
microscopy 

The immunostaining procedure was performed as previ- 
ously described (29). The primary antibodies were as 
follows: Rabbit anti-LHX2 (1:200, AB5756, Chemicon); 
Mouse anti-PAX6 (1:50, Hybridoma bank); Mouse anti- 
NESTIN (1:100, AB5326, Chemicon); Mouse anti-Tujl 
(1:200, MAB1637, Chemicon); Mouse anti-NEUN (1:50, 
MAB377, Chemicon); Mouse anti-OCT4 (1:50, sc-9081, 
Santa Cruz); Mouse anti-SSEA4 (1:200, MAB4304, 
Millipore); Goat anti-OTX2 (1:200, AF1979, R&D), 
and; Goat anti-CERl (1:50, sc-15134, Santa Cruz). The 
secondary antibodies were as follows: Alexa594-donkey 
anti-mouse (Molecular Probe); Alexa488-donkey anti- 
rabbit (Molecular Probe); Alexa594-donkey anti-goat 
(Molecular Probe). For histochemical analysis, H&E 
staining was performed using an Autostainer XL Leica 
ST5010 (Leica Microsystem). The number of rosette-like 
structures per section was counted within 60 out of 200 
serial sections in 1000 mm 3 teratomas, n = 3 (Figure 3F). 
All images were recorded with a Leica TCS SP5 confocal 
microscope. 

Flow cytometry 

Cells were trypsinized, and then resuspended in 0.1% 
paraformaldehyde-phosphate buffered saline (PBS) for 
lOmin on ice. Cells were subsequently transferred to 
90% methanol for 30min on ice. After blocking with 
2% normal donkey serum-PBS, cells were incubated 
with diluted primary antibodies for 1 h at room tempera- 
ture. The primary antibodies were as follows: Rabbit anti- 
LHX2 (1:200, AB5756, Chemicon) and Mouse anti-PAX6 
(1:50, Hybridoma bank). After washing with 0.05% 
Tween 20-PBS, cells were incubated with diluted second- 
ary antibodies for 30min at room temperature. The sec- 
ondary antibodies were as follows: Alexa488-donkey anti- 
mouse (1:500, Molecular Probe) and Alexa488-donkey 
anti-rabbit (1:500, Molecular Probe). After washing 



three times with 0.05% Tween 20-PBS, cells were sus- 
pended in 2% normal goat serum, 0.1% sodium azide 
and 0.05% Tween 20-PBS .The results were recorded 
and analyzed with a BD FACSCalibur flow cytometer. 

Western blot 

Western blotting was performed as previously described 
(29), and separate blots were used for individual 
antibodies. The primary antibodies were as follows: 
Rabbit anti-LHX2 (1:2000, AB5756, Chemicon); Mouse 
anti-PAX6 (1:500, Hybridoma bank); Mouse anti-SOXl 
(1:1000, LS-C39194, Life span); Rabbit anti-SMADl 
(1:1000, 6944, Cell Signaling); Rabbit anti-phosphor- 
SMAD1 (1:1000, 9511, Cell Signaling); Rabbit 
anti-SMAD2 (1:1000, 3102, Cell Signaling); Rabbit anti- 
phosphor-SMAD2 (1:1000, 3101, Cell Signaling); Rabbit 
anti-P-catenin (1:1000, 9581, Cell Signaling); Rabbit anti- 
phosphor- (3-catenin (Ser33/37/Thr41) (1:1000, 9561, Cell 
Signaling); Mouse anti-P-Actin (1:100000, A5441, 
Sigma). The secondary antibodies were as follows: HRP- 
donkey anti mouse (1:6000, Sigma); HRP-goat anti rabbit 
(1:6000, Sigma); HRP-donkey anti goat (1:3000, Jackson). 
Signals were detected and recorded with an ECL detection 
system (Thermo) and a luminescent image analyzer (LAS- 
4000 mini, Fujifilm). 

Microarray analysis 

Total RNA was extracted from iLHX2 cells at day 4 and 
day 16 after doxycycline treatment under conditions un- 
favorable to neural differentiation (3% FBS-DMEM), 
and from shLHX2 cells at day 12 of IVD in N2 media 
(11), using TRIzol reagent (Molecular Research Center). 
Two biological replicates per condition were used. All 
gene expression results were obtained by the Illumina 
Gene Expression Service Laboratory at National Taiwan 
University, Taiwan, using the Illumina microarray 
platform. A HumanHT-12 v4 Expression BeadChip was 
scanned with a BeadStation 450 Scanner, and data were 
analyzed using GeneSpring X software (Agilent) and 
Ingenuity Pathway Analysis (http://www.ingenuity.com/). 
Raw data were normalized independently for each experi- 
ment using Robust Multichip Average. Three pairwise 
comparisons were performed with GeneSpring software: 
day 4 doxycyline-treated and non-treated iLHX2 cells; 
day 16 doxycyline-treated and non-treated iLHX2 cells; 
day 12 differentiating shLuc- and shLHX2-expressing 
hESCs. Changes in gene expression found to have an 
adjusted P < 0.05 and a fold change >2 were considered 
to be significantly different. Gene Ontology (GO) enrich- 
ment was determined by annotating each gene with the 
relevant lineage (ectoderm, mesoderm or endoderm); 
neural-related genes were identified by GO descriptions 
on the NCBI Web site (http://www.ncbi.nlni.nih.gov/ 
gene). For signal pathway analysis, the results of two 
groups (day 4 doxycyline-treated and non-treated iLHX2 
cells, and day 12 differentiating shLuc- and shLHX2-ex- 
pressing hESCs) were analyzed using Ingenuity Pathway 
Analysis Software. Genes were annotated with signal 
pathways related to neural development (Wnt/P-catenin, 
TGFp, FGF, BMP, Notch and Sonic Hedgehog signal 
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pathways) using Ingenuity Pathway Analysis Software. 
Pathways with P < 0.01 [— log(T value) >2] were considered 
significantly different. The NCBI accession number for the 
microarray data reported in this article is GES31810. 

Chromatin immunoprecipitation 

Chromatin immunoprecipitation (ChIP) was performed as 
previously described (31). Briefly, undifferentiated H9 
hESCs and H9 derived neural precursors isolated from 
rosettes at day 16 of IVD were treated with 1% (v/v) 
formaldehyde for 8min at room temperature, and the 
cross-linking reaction was then quenched with 125 mM 
glycine. Chromatin was sonicated to an average size of 
100-150 bp. Each reaction contained 5 x 10 5 cells and 2ug 
antibody (goat anti-LHX2 antibody, SC- 19344, Santa Cruz 
Biotechnology; rabbit anti-H3K4mel antibody, rabbit anti- 
H3K27ac antibody, Abeam). The isolated DNA was 
purified using a QIAquick PCR purification kit (Qiagen). 
The purified DNA was quantified by qPCR using SYBR® 
FAST 2X qPCR Master Mix (KAPA) and an ABI 
PRISM© 7900 Sequence Detection System (Applied Biosy- 
stems), with H9 hESC DNA as standard. Enrichment (% 
Input) was determined as ChIP DNA/ input DNA x 100%. 
Region positions are listed in Supplementary Table S2 and 
primers are listed in Supplementary Table SI. Primers for 
qPCR were confirmed to be specific by the presence of a 
single band on electrophoresis on a 2% agarose gel, and of a 
single peak in the dissociation curve. Each experiment was 
repeated at least three times, with each replicate producing a 
similar pattern. The results are expressed as means ± SE of 
three qPCR triplicates. 

Enhancer reporter assay 

PAX6 P0 (1.5 kb) (32) and PI (1.4 kb) (33) promoters, 
CER1 promoter (0.85 kb) (34) and putative enhancers 
(region of each enhancer: primer —1 to +1) (listed in 
Supplementary Table S2) were cloned into pGL3 
(Promega) to generate reporter plasmids, as shown in 
Figures 4E and 7G. pGal containing the P-galactosidase 
gene was used as an internal standard to control for trans- 
fection efficiency. 

Reporter plasmid (0.5 ug) and pGal (0.1 ug) were 
co-transfected into iLHX2 undifferentiated cells with 
Lipofectamine 2000 (Invitrogen), following the standard 
Lipofectamine 2000 transfection protocol. Two days after 
transfection, the iLHX2 cells were separated from the 
feeder cells, and treated with 2(ig/ml doxycycline. After 
doxycycline induction for 3 days, cells were lysed with 
Reporter Lysis Buffer (Promega). Luciferase activity was 
examined using the Luciferase Assay System (E4030, 
Promega) according to the manufacturer's standard 
protocol, and luminescence was measured using a 
Luminescence Counter (PerkinElmer). The internal 
control ((3-galactosidase) was examined using the 
P-galactosidase Enzyme Assay System (E2000, Promega) 
according to the manufacturer's standard protocol, and 
signals were measured by TopCount NXT Microplate 
Scintillation. Relative luciferase units (RLU) were 
calculated as the activity of luciferase/P-galactosidase 
normalized to cells untreated with doxycycline. 



Statistical analyses 

All experiments were performed in triplicate and the 
results are shown as the mean ± SE. Student's / test was 
used for statistical analysis; P < 0.05 was considered 
significant. 

RESULTS 

LHX2 is one of the earliest markers for hESC-derived 
neural progenitors 

To investigate the expression pattern of LHX2 in human 
neural differentiation, H9 hESCs were plated and cultured 
in neural differentiation media (Figure 1A), and then sub- 
jected to qRT-PCR and western blotting analysis of 
LHX2 and various neural markers. LHX2 mRNA was 
readily induced in hESCs at day 4 of differentiation, and 
its expression preceded the induction of other neural pro- 
genitor marker genes, such as PAX6, SO XI, LM03, 
LIX1, MEIS2 and RAX, and the neuronal gene MAPI 
(Figure IB). Expression of LHX2 plateaued at IVD day 
12, as demonstrated by qRT-PCR (Figure 1C). Western 
blotting analysis revealed that the expression level of 
LHX2 protein was also significantly increased at IVD 
day 4, before the upregulation of PAX6 (Figure ID). To 
examine the spatial expression pattern of LHX2, we 
stained the differentiating hESCs with antibodies against 
LHX2 and other neural markers. Immunocytochemical 
staining revealed that LHX2 was co-expressed with 
PAX6, Nestin or OTX2 in the cells within or around the 
emerging early neural rosette-like structures (Figure IE, 
Supplementary Figure S1A and B). The number of 
LHX2+ and LHX2+/PAX6+ cells increased as differenti- 
ation progressed (Supplementary Figure S1E), and most 
of the PAX6 + cells also expressed LHX2 at IVD day 10 
(Figure IE). After 4 weeks of IVD, certain cells at the 
periphery of the neural rosette structures started to lose 
LHX2 and/or PAX6 expression (Figure IE and 
Supplementary Figure S1A), and began to exhibit the 
morphological features of neurons (Supplementary 
Figure SIC and D). After further differentiation, these 
cells expressed the mature neuronal markers Tujl and 
NeuN; LHX2 was found to be expressed in some of the 
Tujl + (Supplementary Figure SIC) and NeuN + neurons 
(Supplementary Figure SID). This is in agreement with 
previous reports that Lhx2 is expressed in a subpopulation 
of postmitotic neurons (18,35). Together, these results 
demonstrate that LHX2 expression starts at an early 
stage of neural induction in differentiating hESCs, and 
precedes the expression of PAX6 and SOX1 in hESC-NPs. 

LHX2 promotes neural differentiation 

Next, we performed conditional gain-of-function studies 
to investigate if expression of LHX2 is sufficient for neural 
fate induction in hESCs. To this end, we established 
hESCs that ectopically overexpressed LHX2 (iLHX2 
hESCs) by using electroporation and antibiotic selection 
to introduce the inducible Tet-On constructs TetR and 
TetO-LHX2 into undifferentiated H9 hESCs (Supple- 
mentary Figure S2). After doxycycline induction (dox + ) 
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for 4 days (Figure 2A), the differentiating iLHX2 hESCs 
contained a large number of emerging Nestin + cells and 
rosette-like structures containing LHX2 + and PAX6 + , 
whereas iLHX2 hESCs without doxycycline induction 
(dox~) maintained the typical morphological traits of 
hESCs without rosette-like structures (Figure 2B). 
Furthermore, immunocytochemical and flow cytometric 
analyses demonstrated that the number of LHX2- and 
PAX6-expressing cells was significantly increased in the 
dox + iLHX2 hESCs, as compared with the dox~ iLHX2 
hESCs (Figure 2B and C). The increased level of PAX6 
protein in the LHX2-overexpressing cells was also con- 
firmed by western blot analysis (Figure 2D). Furthermore, 
qRT-PCR analysis showed that expression of the 
pluripotency gene, NANOG, was decreased, whereas ex- 
pression of the early neural precursor marker genes PAX6, 
SOX1, SIX6 (36,37), LM03 (38,39), LIX1 (39,40), MEIS2 
(41,42) and NCAD (8,43) [which are highly expressed in 
hESC-derived neural progenitors (7,8) and upregulated by 
PAX6 (11)], and the regional neural genes, including fore- 
brain gene FOXG1, midbrain gene PITX3 and hindbrain 
gene HOXB4, were significantly increased in dox + as 
compared with dox~ iLHX2 hESCs (Figure 2E). Similar 
effects of LHX2 overexpression on neural gene expression 
were also observed in several other hESC lines and iPSC 
lines (Supplementary Figure S3). Together, these data 
indicate that LHX2 overexpression is sufficient to 
promote and facilitates neural differentiation in pluripo- 
tent stem cells. 

Disruption of LHX2 expression compromised neural 
differentiation in hESCs 

To examine whether LHX2 is required for neural differ- 
entiation in hESCs, LHX2 was knocked down by consti- 
tutive overexpression of LHX2 shRNA. Two lentiviral 
shRNA vectors targeting the LHX2 homeodomain 
sequence, or shLuc as a control targeting the luciferase 
gene, were introduced into hESCs, and puromycin selec- 
tion was performed to establish stable hESC clones 
(Supplementary Figure S2). After neural induction for 
16 days, distinct neural rosette-like structures were 
found in shLuc-expressing hESCs, whereas the shLHX2- 
expressing hESCs retained undifferentiated phenotypes 
(Figure 3A). In addition, the expression level of LHX2 
and PAX6 proteins, and the number of LHX2 and 
PAX6-expressing cells, were significantly decreased in 
differentiating shLHX2-expressing hESCs, as revealed by 
western blot and flow cytometry (Figure 3B and C), re- 
spectively. To examine whether knockdown of LHX2 
affected the expression of pluripotency and neural genes, 
qRT-PCR analysis was performed on RNA extracted 
from shLuc and shLHX2-expressing hESCs at day 0 and 
day 12 of neural differentiation (Figure 3D). shLHX2-ex- 
pressing hESCs presented with significantly upregulated 
expression of the pluripotency genes OCT4 and 
NANOG, and significantly reduced expression of neural 
progenitor genes, such as PAX6, and regional neural 
genes, especially, forebrain associated genes, as 
compared with shLuc-expressing hESCs (Figure 3D). 



In addition to the above in vitro experiments, we also 
examined whether knockdown of LHX2 affected neural 
differentiation in vivo by teratoma formation. To generate 
teratomas, hESCs constitutively expressing shLHX2 or 
shLuc were subcutaneously injected into NOD/SCID 
mice. There was no difference in the efficiency of generation 
or the weight of the LHX2 knockdown (shLHX2) or 
control (shLuc) hESC-derived teratomas (data not 
shown). By subjecting a section of teratomas to H&E 
staining, we found that neural rosette-like structures were 
significantly more abundant in teratomas derived from 
shLuc-expressing hESCs than in teratomas derived from 
shLHX2-expressing hESCs (Figure 3E and F). 
Importantly, the levels of LHX2, PAX6 and SOX1, and 
the numbers of LHX2 + , PAX6 + and Nestin + cells in tera- 
tomas derived from shLHX2-expressing hESCs were also 
reduced as compared with shLuc-expressing hESCs (Figure 
3E and G). Taken together, these results suggest that LHX2 
is critical for the neural differentiation of hESCs. 

LHX2 regulates PAX6 expression by binding to its 
enhancers during neural differentiation 

Our gain- and loss-of-function studies revealed that LHX2 
is both necessary and sufficient for neural differentiation 
in hESCs. We therefore went on to investigate the molecu- 
lar mechanisms by which LHX2, in its role as a transcrip- 
tion factor, controls neural differentiation. Based on the 
earlier finding that PAX6 is a human neuroectoderm de- 
terminant (11), and our finding that LHX2 expression 
precedes that of PAX6 (Figure 1), we postulated that 
LHX2 may regulate PAX6 to direct neural differentiation 
in hESCs. To test this hypothesis, we first examined 
whether LHX2 binds to the PAX6 promoters directly. 
To this end, we screened for LHX2-binding sites (TAAT 
TA) (44) in the PAX6 promoters P0 and PI (responsible 
for the expression of PAX6a and PAX6a and b, respect- 
ively). However, we did not find complete LHX2 binding 
motifs in the PAX6 promoters (Figure 4A and 
Supplementary Figure S4). Furthermore, ChIP and 
luciferase reporter assays showed that LHX2 does not 
bind to either of the PAX6 promoters (Supplementary 
Figures S4C) and has no effect on the activities of the 
PAX6 promoters (Figure 4E). 

Because recent studies have shown that enhancers are 
highly cell-type specific (45^17), we further tested whether 
LHX2 binds to PAX6 enhancers (48-52). We screened for 
potential LHX2-binding sites by searching for LIM HD- 
sensitive regulatory regions, which contain a LHX2 con- 
sensus binding sequence (TAATTA) close to a second 
LIM HD-binding sequence (TTAATAAA) (44); this 
resulted in the identification of three PAX6 enhancers 
containing putative LHX2-binding sites (Figure 4B). 
Furthermore, ChIP using an anti-LHX2 antibody and 
chromatin prepared from undifferentiated hESCs and 
hESC-NPs demonstrated that LHX2 bound to sites 
residing in these three PAX6 enhancer regions (E-6k, 
ETel and E+156k) in hESC-NPs (Figure 4C and 
Supplementary Table S2). As it has been recently 
reported that mono-methylation of H3 lysine 4 
(H3K4mel) and acetylation of H3 lysine 27 (H3K27ac) 
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Figure 2. Conditional overexpression of LHX2 promotes the transition of hESCs into the neuroectoderm lineage. (A) Procedure of doxycycline 
induction in iLHX2 hESCs. (B) Immunocytochemical analysis of iLHX2 hESCs treated with (dox + , lower panel) or without (dox - , upper panel) 
doxycycline for 3 days with specific antibodies against LHX2, PAX6 and Nestin. Nuclei were counterstained with DAPI (blue). (C) Flow cytometric 
analysis of a number of LHX2 + or PAX6 + cells during IVD of iLHX2 hESCs treated with doxycycline for 3 days. Gating was determined using H9 
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(right panels, = 4). (D) Western blot of proteins isolated from iLHX2 hESCs after doxycycline induction for 3 days using specific antibodies against 
LHX2 and PAX6. Differentiating H9 ESCs and dox~ iLHX2 hESCs were used as controls. (E) qRT-PCR for the indicated genes using mRNA 
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Figure 3. Reduction of LHX2 expression impaired in vitro and in vivo neural developmental potential of hESCs. (A) Morphology and immuno- 
cytochemical analysis of differentiating shLuc- (upper panel) or shLHX2-expressing (lower panel) hESCs at IVD day 16 using antibodies against 
LHX2 and PAX6. Nuclei were counterstained with DAPI (blue). (B) Western blot analysis of LHX2 and PAX6 protein expression in shLHX2 and 
shLuc-expressing hESCs at the indicated days. Undifferentiated hESCs were used as a negative control. (C) Flow cytometric analysis of the number 
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negative control. The average percentage of LHX2 + or PAX6 + cells is indicated (right panels, n = 3). (D) qRT-PCR for the indicated genes using 
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serve as chromatin signatures of active enhancers and can 
predict the developmental state of a given enhancer (46,47), 
we next used ChlP to investigate the appearance of these 
chromatin marks at LHX2-binding sites in PAX6 enhan- 
cers in hESCs and hESC-NPs (Figure 4D). We found that 
the LHX2-binding sites in E-6kb, ETel (site A) and 
E+156kb (site A and C) enhancers were all enriched with 
both H3K4mel and H3K27ac marks in hESC-NPs, but 
not in hESCs. This suggests that the LHX2 binding sites 
on PAX6 enhancers may be inactive in ESCs and activate 
in neural progenitors, and as such these signatures may be 
indicative of a direct role for LHX2 in activating PAX6 
transcription during neural development. 

To further examine whether these newly identified 
neural-specific enhancers play a direct role in controlling 
PAX6 expression, we set up an enhancer reporter assay in 
inducible LHX2 overexpression hESCs (iLHX2 hESCs). 
The LHX2-binding sites in E-6kb, ETel (site A) and 
E+156kb (site A and C) enhancers (Figure 4C) were 
used for the luciferase reporter assay. The reporter con- 
structs were generated by fusing a luciferase reporter with 
the selected enhancers and their corresponding promoters 
(48,52) (Figure 4E), and were then introduced into iLHX2 
hESCs. Reporters containing the ETel-A and E+156kb-A 
PAX6 enhancer regions resulted in significantly greater 
luciferase activity in dox + iLHX2 cells than reporters con- 
taining the promoter alone (Figure 4E). However, 
mutation of the LHX2-binding motif in these enhancers 
abolished the effect (Figure 4E). Together, these results 
suggest that LHX2 is recruited to active PAX6 enhancers, 
and promotes PAX6 expression during neural lineage 
differentiation. 

LHX2 regulates neural differentiation through PAX6 

The ChlP and enhancer luciferase assays demonstrated 
that LHX2 modulates the neural differentiation of 
hESCs through regulating PAX6 expression; however, 
whether LHX2 primarily targets PAX6 to induce neural 
differentiation in hESCs is unclear. We hypothesized that 
if LHX2 promotes neural differentiation primarily 
through PAX6, a reduction of PAX6 expression in 
hESCs would compromise LHX2-mediated neural differ- 
entiation. To test this hypothesis, we manipulated expres- 
sion of LHX2 and PAX6 in hESCs simultaneously, and 
investigated the effects on morphology and the expression 
profile of a panel of neural genes, including SOX1, LM03, 
SIX6, LIX1 and NCAD (as determined by qRT-PCR). 
PAX6 knockdown overrode the effect of LHX2 
overexpression on both neural rosette formation and 
neural gene induction (Figure 5A). On the contrary, the 
same set of neural genes was upregulated, and abundant 
neural rosettes could still be observed when PAX6 was 
overexpressed in LHX2 knockdown cells (Figure 5B). 



Therefore, these results demonstrate that LHX2 is an es- 
sential neural transcription factor upstream of PAX6, and 
LHX2 primarily targets PAX6 to regulate the neural 
differentiation of hESCs. 

LHX2 affects transcription factor expression and 
signaling pathways involved in neural differentiation 

To determine the effect of manipulating LHX2 expression 
on the global gene expression profile of hESCs, we per- 
formed transcriptome analysis using a cDNA microarray, 
and mRNA extracted from iLHX2 cells treated with or 
without doxycycline to induce LHX2 overexpression. 
Using the fold change analysis tool of GeneSpring 
software, 209 genes at day 4 and 661 genes at day 16 
with at least a 2-fold average expression change between 
dox~ and dox + iLHX2 hESCs were selected and subjected 
to GO analysis. Intriguingly, the majority of the neural 
genes in dox + iLHX2 hESCs showed increased expression 
at both day 4 and day 16 (Figure 6A). Some of the neural 
genes with increased expression at day 16 in dox + iLHX2 
hESCs are postmitotic neuronal genes (MAP2, 
NEUROG2), and others are involved in regional pattern- 
ing of the central nervous system (CNS) (FOXG1, PTH2, 
EMX1, OLIG2). These results confirmed our earlier ob- 
servation in dox+ iLHX2 hESCs that LHX2 induces ex- 
pression of neural genes (Figure 2), and thereafter 
activates the neuronal program. We subsequently per- 
formed transcriptome analysis to investigate the effects 
of LHX2 knockdown in hESCs. We focused on 41 out 
of 86 embryonic germ layer-related genes, which have at 
least a 2-fold average expression change between shLuc- 
and shLHX2-expressing hESCs at day 12 of IVD. 
Interestingly, the expression of most of the neural genes 
was reduced in differentiating shLHX2-expressing hESCs 
(Figure 6A), confirming our observation in Figure 3. 
Overall, the findings of the genome-wide gene expression 
analysis of iLHX2 and shLHX2 hESCs suggest that 
LHX2 expression has a profound impact on transcrip- 
tional regulation during neural differentiation. 

We also investigated whether LHX2 directs neural dif- 
ferentiation by affecting relevant extrinsic signaling 
pathways, in addition to its role in regulating transcrip- 
tional factors. To this end, we performed Ingenuity 
Pathways Analysis on the gene expression profiles of 
hESCs with LHX2 overexpression or knockdown. 
Through this, we identified four pathways affected in 
hESCs by altered LHX2 expression: the WNT/P-catenin, 
TGFp, FGF and BMP signaling pathways (Figure 6B). 
We confirmed downregulation of the WNT/P-catenin, 
TGFp and BMP signaling pathways by western blot: 
phosphor-SMADl (a downstream effector of BMP) and 
phosphor-SMAD2 (a downstream effector of TGFp) were 
downregulated, and phosphor- P-catenin (a downstream 



Figure 3. Continued 

and Nestin (arrows indicate neural rosette-like structures). Nuclei were counterstained with DAPI (blue). (F) Quantitative analysis of the number of 
neural rosette-like structures in teratomas derived from shLHX2 (831 and 968) or shLuc-expressing hESCs (n = 3). (G) Western blot analysis of the 
teratomas derived from shLHX2- or shLuc-expressing hESCs with specific antibodies against LHX2, PAX6 and SOX1. Error bars represent the 
mean ± SEM. Significance: *P<0.05. ShRNA were named according to the targeted gene and the first shRNA targeting nucleotide in the gene 
coding sequence (number). For the targeting sequences of shLHX2-831 shLHX2-968, please see the 'Materials and Methods' section. 
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Figure 4. LHX2 regulates PAX6 expression through its enhancers in hESC-NPs. (A) Schematic representing the genomic organization of the PAX6 
promoter regions, P0 and PI, from which full-length PAX6 mRNA is transcribed. PO spans positions —1450 to —1 and PI spans positions +6630 to 
+8083 with respect to the transcription starting site (TSS). The filled triangles indicate incomplete LHX2 binding sites. ChIP primers used in the 
study are indicated by black dots. '0' indicate regions containing incomplete LHX2 binding sites and T indicates regions located at the 5' end of '0'. 
(B) A schematic representation of the putative enhancers of PAX6. Unfilled boxes indicate PAX6 enhancers containing putative LHX2 binding sites 
(— 6K, ETel and +156K). Filled boxes indicate PAX6 enhancers without putative binding sites (EP, EE, — 2K, E1E and Ele4H). Unfilled triangles 
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Figure 5. PAX6 is a major downstream target of LHX2 for regulation 
of neural differentiation. (A) qRT-PCR analysis of the indicated genes 
using mRNA isolated from hESCs constitutively expressing shLuc or 
shPAX6 and transiently expressing GFP or LHX2 by lentiviral infec- 
tion for 3 days. Bright field images show the morphological change in 
hESCs 6 days after viral infection (compare a and b). The number of 
neural rosette-like structures and the expression level of the indicated 
neural genes were dramatically reduced by PAX6 knockdown. 
(B) qRT-PCR analysis of the indicated genes using mRNA isolated 
from hESCs constitutively expressing shLuc or shLHX2 and transiently 
expressing GFP or PAX6 by lentiviral infection for 3 days (compare a 
and b). Bright field images show the morphological changes in the 
hESCs 6 days after viral infection. The number of neural rosette-like 
structures and the expression level of the indicated neural genes were 
not affected by knockdown of LHX2 in hESCs overexpressing PAX6. 
Error bars represent the mean ± SEM (n = 3). Significance: *P < 0.05. 
Dotted lines indicate rosette-like structures. 



effector of WNT) was upregulated, when LHX2 was 
overexpressed in differentiating hESCs (Figure 6C). We 
then examined the gene expression profiles of hESCs 
with LHX2 overexpression or knockdown to identify the 
genes through which LHX2 affects these signaling 
pathways (Figure 6D). Based on GO terms, extracellular 
components and transmembrane activity, we identified 
four candidate genes with expression strongly correlated 
to that of LHX2: CER1 (encoding a BMP, nodal and 
WNT antagonist), TWSG1 (encoding a BMP antagonist), 
SFRP2 (encoding a negative regulator of the WNT recep- 
tor) and LEFTY 1 (encoding a TGFp antagonist) 
(Figures 6D and 7A). 

LHX2 attenuates BMP and WNT signaling and inhibits 
non-neural differentiation by augmenting expression of 
CER1 

To test our hypothesis that LHX2 promotes neural differ- 
entiation by suppressing the non-neural program through 
targeting signaling antagonists, we examined whether 
LHX2 directly targets CER1, the expression of which 
was most significantly affected by manipulation of 
LHX2 expression (Figure 7A). We first examined the 
temporal expression pattern of CER1 during neural dif- 
ferentiation using RT-PCR and qRT-PCR. We found that 
CER1 expression was initiated in differentiating hESCs at 
day 6 of IVD, and peaked at days 8-10, when neural 
progenitors were evidently emerging in the culture 
(Supplementary Figure S5A); this suggests that expression 
of LHX2 is followed by that of CER1. In addition, CER1 
was observed to be co-expressed with LHX2 at day 10 and 
day 20 of IVD, suggesting that initiation of CER1 expres- 
sion is highly correlated with that of LHX2 during neural 
differentiation of hESCs (Supplementary Figure S5B). 

To establish whether CER1 plays a role in human 
neural differentiation, loss-of-function studies were con- 
ducted with hESCs constitutively expressing shCERl. 
The use of qRT-PCR revealed that CER1 knockdown 
significantly decreased expression of both neural progeni- 
tor genes (SOX1, SIX6, LM03, LIX1, MEIS2 and 
NCAD) and other CNS regional-specific genes [with the 
exception of midbrain-specific genes (PITX3 and 
LMXla), a hindbrain-specific gene (KROX20) and a 
spinal cord-specific gene (LHX3)] (Figure 7B); these 
results suggest that knockdown of CER1 in hESCs 
affects the differentiation of multiple neural lineages, 
although its role may be more pronounced in forebrain 
differentiation. In addition, mesoderm- and endoderm- 
related genes (MIXL1, T, GSC and SOX 17) were 
increased at day 12 of IVD in shCERl -expressing 
hESCs, suggesting that CER1 may function as a repressor 



Figure 4. Continued 

regions located at the 3' end of '0'. (C) Combined ChIP and qPCR analysis of DNA isolated from hESC-NPs using a specific antibody against 
LHX2 revealed that LHX2 exhibits enhanced binding to certain PAX6 enhancers in hESC-NPs, as compared with undifferentiated hESCs (DNA 
from the latter was used as a control). (D) Combined ChIP and qPCR analyses of hESC-NP DNA using specific antibodies against H3K4mel (upper 
panel) or H3K27ac (lower panel) revealed enhanced occupancy of these chromatin marks on LHX2 binding sites of certain PAX6 enhancers. 
(E) Mutational analysis of LHX2 binding sites in PAX6 enhancers in iLHX2 hESCs with or without doxycycline induction. Reporter constructs 
containing wild-type or mutated LHX2 binding sites in enhancer —6k (E-6k), ETel-A (ETel-A) or +156K (E+156K-A and -C) and PAX6 promoters 
P0 and PI were generated, and used to evaluate the role of the enhancers in regulating activation of PAX6 promoters and PAX6 expression under 
LHX2 overexpression (LHX2+). Error bars represent the mean ± SEM (n = 3). Significance: *P<0.05. 
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Figure 6. Transcriptome profiling of hESCs with altered LHX2 expression revealed that LHX2 affects the expression of neural genes and signaling 
pathways. (A) The heat map shows filtered neural genes with expression significantly affected by manipulation of LHX2 expression in hESCs. 

(B) Identification of signaling pathways affected by manipulation of LHX2 expression in hESCs by Ingenuity Pathway Analysis software. P<0.01 
was used as the cutoff value for selecting neural signaling pathways significantly affected by manipulation of LHX2 expression in hESCs. 

(C) Western blot analysis of proteins isolated from dox + iLHX2 hESCs with antibodies against ACTB and proteins downstream of the BMP4 
(phosphor-SMADl and SMAD1), TGFp (phosphor-SMAD2 and SMAD2) and WNT (phosphor- P-catenin and P-catenin) signaling pathways. 
Proteins isolated from dox~ iLHX2 hESCs were used as a control. (D) The heat map shows genes annotated to the indicated signaling pathways 
by Ingenuity Pathway Analysis software with increased expression in dox + iLHX2 cells (fold change >1.3), and decreased expression in 
differentiating shLHX2-expressing hESCs (fold change >1.3). Green and red values represent fold change for down- and upregulation, respectively 
(A and D). The fold change of gene expression was calculated using GeneSpring software as iLHX2 hESC with doxycycline induction to without 
(iLHX2 hESC dox + /dox~) or shLHX2- to shLuc-expression hESCs (shLHX2/shLuc) on the days indicated. 
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Figure 7. LHX2 attenuates BMP and WNT signaling by augmenting CER1 expression. (A) qRT-PCR analysis of the expression of CER1, TWSG1, 
SFRP2 and LEFTY1 genes using mRNA isolated from iLHX2 hESCs (iLHX2-3 and -18, induced by doxycycline for 3 days) or shLuc- or shLHX2- 
expressing hESCs (IVD day 12). ACTIN was used as the internal control for normalization. Fold change was normalized independently using iGFP 
expression or shLuc expression in parallel as a control (n = 3). (B) qRT-PCR analysis of indicated genes using mRNA isolated from differentiating 
shCERl-expressing hESCs at IVD day 12 (n = 3). (C) Immunocytochemical analysis using specific antibodies against CER1 and GFP in hESCs 
transiently expressing LHX2-eGFP for 3 days by lentivirus infection. CER1 and LHX2 (eGFP) were co-expressed in hESCs. (D) Schematic rep- 
resentation of the identified LHX2 binding sites (arrows, Supplemental Figure S6B) neighboring the CER1 gene. (E) Combined ChIP and qPCR 



(continued) 



7766 Nucleic Acids Research, 2013, Vol. 41, No. 16 



of mesoderm- and endoderm-related genes during IVD of 
hESCs (Supplementary Figure S6A). To confirm that 
CER1 is sufficient to mediate LHX2-dependent regulation 
of both WNT and BMP (Figure 6C), we overexpressed 
CER1 in shLHX2-expressing hESCs, and investigated 
the effect on the downstream proteins of BMP and 
WNT signaling by western blot (Supplementary Figure 
S6B). This revealed that overexpression of CER1 is suffi- 
cient to attenuate BMP and WNT signaling. 

Next, we examined whether overexpression of LHX2 
(through transient infection with lentivirus encoding 
LHX2-eGFP) results in upregulation of CER1 expression 
in hESCs. Immunostaining revealed that CER1 was ex- 
pressed in LHX2-overexpressing cells 4 days after viral in- 
fection, suggesting that LHX2 may regulate CER1 
expression (Figure 7C). These results, combined with our 
finding that PAX6 is the primary downstream target of 
LHX2, raise the possibility that LHX2 may regulate 
CER1 expression through PAX6. To address this, we tran- 
siently overexpressed shPAX6 in doxycycline-treated 
iLHX2 hESCs by lentiviral infection (Supplementary 
Figure S5C). CER1 expression continued to be significantly 
increased in iLHX2 hESCs, even on PAX6 knockdown. 
This result suggests that LHX2 may regulate CER1 
through a PAX6-independent pathway in hESCs. 

To elucidate the direct transcriptional mechanisms that 
govern CER1 expression by LHX2, we sought to identify 
cis-acting regulatory elements by aligning the CER1 
region (Chr9:14673331-14769701) of the human genome 
to various species, including rhesus monkey, dog, horse, 
mouse, rat and chicken. Through using the VISTA 
browser to identify the degree of evolutionary constraint 
in the CER1 gene, we identified nine conserved regions, 
seven of which contained putative LHX2 binding sites 
(Supplementary Figure S7B). LHX2-ChIP analysis 
further revealed that LHX2 specifically targeted 5 out of 
20 predicted sites in the 7 conserved regions in hESC-NPs 
(Figure 7D and E), but not the CER1 promoter (34) 
(Supplementary Figure S7A). Furthermore, these regions 
were also enriched with the active enhancer marks 
H3K4mel and H3K27ac in hESC-NPs (Figure 7F). To 
further test whether these five newly identified regions 
serve as putative enhancers to regulate CER1 promoter 
activities, we performed enhancer reporter assays in 
iLHX2 hESCs. We report that three out of the five 
putative enhancers of CER1 were activated by LHX2, 
and acted as enhancers to regulate the CER1 promoter, 
while mutation of the LHX2-binding site eliminated this 
effect (Figure 7G). Collectively, these results suggest that 
LHX2 may also contribute to neural differentiation via 
upregulation of CER1, which itself encodes an antagonist 
of the BMP and WNT/P-catenin signaling pathways. 



DISCUSSION 

Previous studies have reported that the extrinsic signals 
WNT, SHH, BMP and FGF, and the intrinsic factor 
PAX6, regulate neural differentiation (2,3,11). However, 
the transcriptional regulation underlying early neural 
lineage differentiation remains poorly understood. The 
aim of the current study was to increase our understand- 
ing of the transcriptional regulatory network governing 
human early neural lineage differentiation. We provide 
evidence that LHX2 is expressed in the early neural 
lineage derived from human pluripotent stem cells. By 
combining global gene expression profiling with biochem- 
ical analysis of hESCs with altered expression levels of 
LHX2, we found that LHX2 also affects extracellular sig- 
naling pathways, including the WNT, TGF(3 and BMP 
pathways. This reveals that LHX2 is likely to be 
involved in regulating the neural differentiation of 
hESCs through multiple mechanisms. We demonstrate 
here that expression of LHX2 leads to neural fate acqui- 
sition at two levels (Figure 8). First, we report that LHX2 
regulates PAX6 (a transcription factor critical for human 
early neural differentiation) through binding to its enhan- 
cers (enriched with active chromatin markers), and 
through this mechanism LHX2 drives neural lineage dif- 
ferentiation. Second, we have demonstrated that LHX2 
suppresses the WNT and BMP signaling pathways in 
differentiating hESCs through regulation of CER1 expres- 
sion; CER1 acts as an antagonist to attenuate the BMP 
and WNT signaling pathways. Therefore, we have 
identified LHX2 as a key factor for regulating the inter- 
play between intrinsic and extrinsic mechanisms required 
for neural differentiation of hESCs. 

We demonstrate here that LHX2 binds to enhancers, 
rather than promoters, of both PAX6 and CER1. All of 
the identified LHX2 binding sites in the PAX6 (E-6k, 
ETel-A, E+156k-A and C) and CER1 enhancers (A, B, 
C, D and E) are marked by H3K4mel and H3K27ac, 
chromatin marks of active enhancers. These findings are 
in agreement with a previous study that indicated that 
human enhancers are enriched with active chromatin sig- 
natures, and that enhancers, not promoters, reflect cell- 
type-specific gene expression (45). We also noted that 
two of the identified LHX2 binding sites in PAX6 enhan- 
cers, E-6k and E+156k-C, are not active in hESC-NPs. It 
has been previously reported that multiple PAX6 pro- 
moters are regulated by diverse enhancers to transcribe 
different PAX6 isoforms in a cell-type-specific manner. 
For example, the enhancer located 6kb upstream from 
the Pax6 P0 promoter (enhancer —6k) is conserved in 
mice and humans, and regulates P0 promoter activity 
for transcription of PAX6a mRNA in neuroretina (49); 



Figure 7. Continued 

analyses of DNA isolated from hESC-NPs and hESCs with a specific antibody against LHX2 revealed that LHX2 binds to the indicated sites 
neighboring the CER1 gene in hESC-NPs (n = 3). (F) Combined ChIP and qPCR analyses revealed enhanced occupancy of H3K4mel (upper panel) 
and H3K27ac (lower panel) at LHX2 binding sites in the indicated CER1 enhancers in hESCs-NPs, as compared with hESCs (« = 3). (G) Reporter 
constructs containing the indicated wild-type or mutant LHX2 binding sites and CER1 promoters were generated, and used to evaluate the role of 
the enhancers in regulating CER1 promoter activation and CER1 gene expression in iLHX2 hESCs (n = 3). '0' indicates regions containing complete 
LHX2 binding sites. T and '2' indicate regions located at the 5' end of '0' and '—2' and '— V indicate regions located at the 3' end of '0'. Error bars 
represent the mean ± SEM. Significance: */ > <0.05. 
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Figure 8. Model for LHX2-mediated neural differentiation in hESCs. 
LHX2 directly regulates PAX6 expression by binding to its enhancers, 
and promotes neural differentiation of hESCs. On the other hand, 
LHX2 can also regulate CER1 expression by directly binding to the 
putative enhancer region of CER1 and attenuating Wnt and BMP sig- 
naling, further suppressing the mesodermal and endodermal potential 
of hESCs. 



an enhancer located in intron 1 (ETel), which is conserved 
in fugu and humans, regulates PI promoter activity for 
transcription of PAX6a and PAX6b mRNA in telenceph- 
alon, lens and pancreas (48); and enhancers located 156kb 
distal from PAX6 PO promoter (enhancer +1 56k), which 
are conserved in mice and humans, regulate PI promoter 
activity for transcription of PAX6a and PAX6b mRNA in 
the central nervous system (52). Therefore, the likely 
function of the ETel- A and E+156k-A enhancers is to 
regulate neural differentiation in hESCs, and the role of 
the other two identified LHX2-binding enhancers (E-6k 
and E+156k-C) may involve later regional specification. 

Several studies have demonstrated the important role of 
Lhx2 in regulating the extracellular signaling pathways 
involved in neural development. For example, Lhx2 has 
been suggested to couple with BMP signaling to promote 
optic cup ventral identity in mice (21). Recently, it has 
been suggested that Lhx2 may play a role in the regional- 
ization of thalamus formation, through attenuating WNT 
signaling emanating from the epithalamus in zebrafish 
(53). Additionally, expression of Lhx2 in the presumptive 
hippocampus during mouse neocortex formation attenu- 
ates WNT and BMP signaling emanating from the cortical 
hem, thereby restricting cortical hem extension and there- 
after forming a sharp border between the hippocampus 
and cortical hem (18). These observations establish Lhx2 
as an instructive factor that attenuates various signaling 
pathways during forebrain formation. However, the 
mechanism through which Lhx2 regulates its associated 
signaling pathways to control forebrain formation is not 
clear. We identified and confirmed that expression of 



CER1, a gene encoding a secretory ligand that binds to 
WNT and BMP to inhibit the downstream signaling 
cascade, is directly regulated by LHX2 in differentiating 
hESCs. These findings therefore provide a possible mech- 
anism by which Lhx2 regulates downstream signaling and 
controls forebrain patterning. 

The dual roles of Cerl in neural induction and suppres- 
sion of mesodermal or endodermal lineages are well docu- 
mented (54-57). In mice, Cerl is expressed in the anterior 
mesoendoderm, and antagonizes the BMP signal emitted 
from the ventrolateral mesoderm to promote the onset of 
neurulation at early gastrulation (57). Additionally, Cerl 
was also found to induce expression of anterior neural 
genes, such as Ncam and Otx2 in frogs (57). Moreover, 
ectopic expression of CER1-S (a truncated form of CER1, 
which can only target nodal proteins) in hESCs resulted in 
early onset of neural rosette formation, and an increase in 
neural gene expression at the expense of non-neural 
lineages, such as the mesoderm and endoderm (56). In 
line with these observations, we found that several 
mesoderm- and endoderm-related genes are increased in 
differentiating CER1 knockdown hESCs (Supplementary 
Figure S6A), and WNT and BMP signaling is attenuated 
(Supplementary Figure S6B) and neural genes are 
activated in CERl-overexpressing hESCs (data not 
shown). These results indicate that upregulation of 
CER1 by LHX2 may be responsible for the activation of 
the neural program and suppression of the non-neural 
program during early hESC differentiation. 

In addition to CER1, our results indicate that expres- 
sion of TWSG1 (encoding a BMP antagonist), SFRP2 
(encoding a WNT antagonist) and LEFTY 1 (encoding a 
TGFp antagonist) are affected by LHX2 expression levels 
in hESCs. WNT, BMP and TGFfS signaling have been 
previously implicated in the repression of neural differen- 
tiation, and inactivation of these pathways is essential for 
neural induction (3). Thus, an analogous mechanism, 
whereby LHX2 activates antagonists of WNT and BMP 
proteins in certain areas to restrict these signaling 
pathways, may explain the earlier observation that the 
presence of LHX2 in the developing forebrain prevents 
excessive expansion of misplaced tissue (to enable appro- 
priate regionalization) (18). 

Through gene-expression profiling and manipulation of 
LHX2 expression in hESCs, we found that expression of 
several neural transcription factor genes, including MEIS2 
and RAX, overlaps with expression of LHX2 and/or 
PAX6 at differentiation stages, and is strongly correlated 
with LHX2 expression. This suggests that MEIS2 and 
RAX are likely to be involved in the LHX2 and/or 
PAX6 transcriptional networks that determine neural 
fate. MEIS2, a member of the TALE family, is expressed 
in the proliferation zone of the human fetal brain, and is 
co-expressed with PAX6 during the development of the 
ganglionic eminence (42). Regulation of PAX6 also con- 
tributes to mouse lens morphogenesis (58). In our gain-of- 
function studies, ectopic LHX2 or PAX6 expression in 
hESCs activated MEIS2 expression (Figure 2), and 
ectopic MEIS2 expression in hESCs increased the expres- 
sion of LHX2 and PAX6 (Hou et al., unpublished obser- 
vation). Thus, it is tempting to suggest that MEIS2 may 
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initially be activated by LHX2 and PAX6, and its expres- 
sion may subsequently serve as a positive feedback cue to 
enhance LHX2 and PAX6 expression. In mice, RAX is 
expressed in the anterior neural plate of E8.5 embryos, 
and is later expressed in the optic cup and ventral fore- 
brain of E10.5 mouse embryos; RAX activates PAX6 ex- 
pression during retina formation (59), suggesting it may 
also play a role in regulating PAX6 expression in neural 
development in mice. Therefore, we propose that RAX 
may also be involved in the LHX2-PAX6-mediated regu- 
latory network controlling neural differentiation. 
Understanding of the functional roles of these factors in 
neural induction/differentiation, and their relationship to 
LHX2 and/or PAX6, will be critical for further dissection 
of the genetic regulatory network controlling neural dif- 
ferentiation in hESCs. 

Although a reduction in LHX2 expression in hESCs 
significantly compromises neural differentiation, we still 
observed expression of neural markers and the formation 
of neural rosette-like structures during differentiation of 
LHX2-knockout hESCs, both in vitro and in vivo 
(teratoma formation). This may be because LHX2 is not 
absolutely required for neural differentiation of hESCs, or 
because shLHX2 does not completely eliminate LHX2 
protein expression (Figure 3B and G). As regards the 
first possibility, (an)other related gene(s) may result in 
some degree of functional redundancy for LHX2 in the 
control of neural differentiation. Such functional redun- 
dancy has been described between Lhx2 and Lhx9 in 
spinal cord development in mice (60). Another possible 
explanation is that LHX2 may only be critical for differ- 
entiation of certain neural cell types, such as cells with 
anterior neural cell identity that constitute the majority 
of IVD neural populations under the current neural induc- 
tion protocol (Supplementary Figure S8A) (7). Our obser- 
vations on forebrain-related genes on the effect of LHX2 
loss of function (Supplementary Figure S8B) further 
support this possibility. Dissection of the distinct role of 
LHX2 in regional neural specification in humans will 
require further examination of its effect on the generation 
of specific neural subtypes from hESCs. 

The kinetics of PAX6 and SOX1 expression during 
neural differentiation are inconsistent across the literature 
(11,61). Such discrepancies may arise from two 
possibilities. The first is the neural induction method, as 
there are two procedures for in vitro neural differentiation 
from hESCs. The key difference between these procedures 
is the presence or absence of SMAD inhibitors, such as 
SB431542 for the TGFp pathway, and Noggin for the 
BMP pathway. Studies that did not use such inhibitors 
(such as the present study) have reported that PAX6 
precedes that of SOX1 (11,39,62). On the contrary, the 
sequential expression of SOX1 and PAX6 was reported 
to be reversed in the presence of inhibitors (61). Further, 
it has been reported that initiation of SOX1 expression 
requires inhibition of the TGFp pathway, and initiation 
of PAX6 expression requires the absence of FGF signaling 
(63). These findings suggest that the kinetic expression of 
PAX6 and SOX1 is influenced by the culture condition. 

Taken together, our findings indicate that LHX2 
mediates neural differentiation not only through intrinsic 



transcriptional regulation of PAX6, but also through sup- 
pression of extrinsic factors of key signaling pathways that 
are essential for early lineage differentiation in human. 

SUPPLEMENTARY DATA 
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Supplementary Tables 1 and 2 and Supplementary 
Figures 1-8. 
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